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A model of 3+1 dimensional gravity with negative cosmological constant coupled 
to abelian gauge fields has been proposed as a gravity dual for Lifshitz like critical 
phenomena in 2+1 dimensions. The finite temperature behavior is described by 
black holes that are asymptotic to the Lifshitz fixed point geometry. There is a 
one-parameter family of charged black holes, where the magnitude of the charge is 
uniquely determined by the black hole area. These black holes are thermodynamically 
stable and become extremal in the limit of vanishing size. The theory also has a 
discrete spectrum of localized objects described by non-singular spacetime geometries. 
The finite temperature behavior of Wilson loops is reminiscent of strongly coupled 
gauge theories in 3+1 dimensions, including screening at large distances. 
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1 Introduction 



AdS/CFT duality [H EJ [3] provides unique access to the physics of supersymmetric 
gauge theories at strong coupling. The original conjecture, and various refinements of 
it, have been tested to the extent that its validity is generally accepted in the theoreti- 
cal high-energy physics community. By now AdS / CFT duality has been generalized in 
many directions and the notion of a gauge theory /gravity correspondence has become 
a standard item in the toolbox of high-energy theory. Recently, considerable effort 
has been put into extending AdS/CFT beyond high-energy physics by constructing 
gravity models that are conjectured to be dual to various condensed matter systems 
[I] -[IS]- In this case there is no supersymmetry, or truncation of higher-dimensional 
supergravity, to justify the duality but this is not necessarily a problem. The gravity 
dual provides a phenomenological description of whatever strongly coupled physics 
that is being modeled, and, as such, it can be useful even if the connection to the 
underlying dynamics cannot be spelled out in detail. Ultimately, the usefulness of a 
dual description of a real physical system is to be judged by its success in explaining, 
or better yet, predicting experimental results. 

In the present paper we develop further a recently proposed gravitational dual 
description [16] of a class of critical phenomena exhibiting unconventional scaling of 
the form 

t -> \% x -> Ax, (1) 
with z 1. The so called Lifshitz theory, 

L = j d 2 xdt ((d t <pf - K (V 2 0) 2 ) ■ (2) 

provides a simple example of a 2+1 dimensional field theory which is invariant under 
precisely this kind of scaling (with z = 2). This, and other related models, have 
been used to model quantum critical behavior in strongly correlated electron systems 

In [TB] it was conjectured that strongly coupled systems with Lifshitz scaling can 
be modeled by a gravity theory with a spacetime metric of the following form 

ds 2 = L 2 (-r 2z dt 2 + r 2 d 2 * + ^ . (3) 

This metric is invariant under the transformation 

t^\% r^^, x^Ax. (4) 
A 

The coordinates (t, r, x 1 , x 2 ) are dimensionless and the only characteristic length scale 
of the geometry is L. It was shown in [TB] that an action coupling four-dimensional 
gravity, with a negative cosmological constant, to a simple complement of abelian 
gauge fields, 

S = J d A x^g (R-2K)- l -J (F m A *F (2) + H {3) A *# (3) ) - c J B {2) A F (2) , (5) 
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can support such a metric with z > 1. Here F( 2 ) = dA^ and ff( 3 ) = dB( 2 ) are a two- 
form and a three-form field strength respectively, and the length scale L is related to 
the cosmological constant A = — 5/L 2 . 

Here we focus on the z = 2 case but we expect that most of our results can be 
generalized to other z values. In section 2, we extend the analysis of [16] to global 
coordinates, and demonstrate the existence of a discrete set of solutions that we call 
Lifshitz stars. We also introduce finite temperature by having a black hole at the 
center of an asymptotically Lifshitz space time. Such a black hole carries an electric 
charge that couples to the two-form gauge field strength. Interestingly, the magnitude 
of the charge is uniquely fixed by the black hole size and these black holes become 
extremal in the small size limit. In section 3 we study the thermodynamic properties 
of black holes in asymptotically Lifshitz spacetime, and conclude that, contrary to 
black holes in asymptotically AdS spacetime, they are thermodynamically stable for 
all black hole sizes. In section 4 we speculate on the nature of the boundary theory, 
and perform a calculation analogous to the evaluation of a Wilson loop in AdS/CFT. 

A number of recent papers have analyzed black hole geometries in gravity duals of 
non-relativistic quantum systems [2T]-[25] but those gravitational models are different 
from the one we consider here, leading to a different spectrum of black holes and 
different thermodynamic properties. 



2 Asymptotically Lifshitz spacetime 

In this section we look for global metrics that solve the equations of motion of (jsl) 
and approach the Lifshitz geometry ^ in an asymptotic limit. We find two types 
of spherically symmetric, static solutions. One is a black hole with a non-degenerate 
event horizon and the other is a smooth geometry that describes a non-singular, 
spherically symmetric concentration of the gauge fields. 

The equations of motion of ^ are easily obtained. The gauge fields satisfy a pair 
of coupled equations, 

d* F {2) = -cH {3) , (6) 
d*H (3) = cF {2) , (7) 

and the Einstein equations are 

Gfj,v - jj2~9^u = -jiF^F* - -g^FxaF^) + -(H^xaH Xa - -g^Hx ap H Xap ). (8) 

The Lifshitz metric (|3| is a solution of the equations of motion if the topological 
coupling between the gauge fields is tuned to be 




We will assume this value for c in what follows. 
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Our ansatz for the metric and the two- and three-form fluxes generalizes the one 
employed in [16] to global coordinates, 



g(r)' 



df 



(10) 



'(2) 



H, 



(3) 



h{r) r A t , 
j{r) r A0q A 00, 
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(16) 



The two- form field strength F^), given by h(r), is an electric field directed radially 
outwards from the origin at r = 0. The topological coupling between Fm and Hc$\ in 
the action implies that the three-form flux is electrically charged and acts as a source 
of the electric field. The above ansatz for thus corresponds to a charged fluid 
whose density is governed by j (r). 

The Einstein equations and the field equations for the gauge fields reduce to a 
system of non-linear first order differential equations, 
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(17) 

(18) 

(19) 
(20) 



In the r>l limit, we can replace the sphere metric, d 2 Q, in (10) by a flat metric, or, 
equivalently, neglect the 1/r 2 terms in equations (17) to (19). In this case we recover 
the Lifshitz solution, / (r) 



g (r) = h (r) = j (r) = 1, as constructed by 
also easy to see that the AdS-Schwarzschild geometry, 
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It is 



(21) 



is a solution of the equations with h(r) = j(r) = but we have not found explicit 
analytic solutions with non-trivial gauge fields. It is, however, straightforward to 
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integrate the system of equations numerically and one can learn a lot about the 
behavior of solutions by a combination of numerics and asymptotic analysis. 



It is convenient to first consider equations (18) - (20), which only involve the three 
functions g(r), h(r), and j(r) and then, given a solution to this system, solve equation 
(17) for the remaining function f(r). We are primarily interested in geometries that 

In order to 



1 as r 



are asymptotically Lifshitz in the sense that f,g,h,j 
study the asymptotic behavior at large r we linearize the system (18) 
the fixed point, 
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(22) 



where g — 1 + Sg, etc. By standard manipulations, the matrix that appears in the 
linear system can be brought into Jordan form, 



-4 






(23) 



from which we read off that at the linear level there are two decaying modes, that 
behave as 1/r 4 and log(r)/r 4 , and a zero mode that does not depend on r. In addition 
to the eigenmodes, the general solution of the linear system includes a universal 1/r 2 



mode that comes from the inhomogeneous term on the right hand side of (22). It is 
straightforward to find the general solution of the linear system and confirm that it 
has all these features, but, since we do not need the details of the full solution for 
what follows, we do not write it down here. 

At large values of r the non-zero modes, including the universal 1/r 2 mode, have 
decayed away leaving only the zero mode behind. The solution then has the form 



S«l + 7, h « 1 + 27, j « 1 + 7, 



(24) 



with 7 <C 1, but this is not the whole story. When non-linear corrections are included, 
the zero mode is lifted and becomes either a marginally growing or marginally decay- 
ing mode. 



To study the evolution of the zero mode at large r, we insert (24) into the original 



non-linear equations of motion (18) - (20), except that, since we are assuming that 
the non-zero modes have already decayed, we drop the 1/r 2 terms in the equations. 
Working to leading non-vanishing order in 7 then gives 



rg (rj 



77 2 , rh'{r) » 27 2 , 



rj (r) 



1 



(25) 



The leading order terms on the right hand side all come with a positive sign. This 
immediately implies that a zero mode of positive amplitude 7 > is a growing mode 
and a solution of the equations of motion where such a mode is present at large r 
cannot be asymptotic to the Lifshitz fixed point. 
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A zero mode of negative amplitude will, on the other hand, slowly decay. The 
decay can be described analytically in a simple fashion. Writing 7(7-) = — 1/logr it 
is straightforward to check that 

g(r) « 1 +7(r) +7 2 (r) + . . . , 

h(r) ss l + 2 7 (r) -7 2 (r) + ..., (26) 
j(r) « 1 + 7(r) - 27 2 (r) + . . . , 

solves the non-linear equations of motion up to terms that are small compared to 
l/(logr) 2 . The solution thus appears to be asymptotic to the Lifshitz fixed point, at 
least as far as the functions g, h, and j go, but the approach to the fixed point is 
extremely slow. In order to decide whether the geometry is truly asymptotic to the 
Lifshitz geometry we have to consider the evolution of the remaining function f(r), 



which enters into the gu component of the metric (10). Inserting a zero mode (24) 



with 7(7*) = —1/logr into the f(r) equation of motion (17) gives 



rf'(r) = -^- + .... (27) 
logr 

This integrates to 

f(r) = (28) 
(logr) 4 

which goes to zero in the r — > 00 limit. A solution with a marginally decaying zero 
mode is therefore not asympotic to the Lifshitz fixed point after all. 

The upshot of all this is that the only solutions of our system that are asymptotic 
to the Lifshitz fixed point are those for which the amplitude of the zero mode of 
the linearized system happens to vanish in the asymptotic region. This requires fine- 
tuning of initial values when solving the non-linear equations of motion, which in turn 
reduces the number of free parameters in the solutions that are of most interest to 
us. 



2.1 Black holes 

We now describe results obtained by integrating the equations of motion numerically. 
We find a one-parameter family of black hole solutions with a non-degenerate horizon 
which are asymptotic to the Lifshitz geometry pj). The characteristic parameter of 
the black hole can be either taken as r , the value of the area coordinate r at the 
horizon, or h = h(r ), the value of the radial electric field at the horizon. At first 
sight, one would expect these two parameters to be independent as they are for an 
ordinary Reissner- Nordstrom black hole but, as discussed above, only a restricted set 
of geometries approaches the Lifshitz fixed point at r —>■ 00. The restriction on the 
parameters can be understood in terms of the interaction between the charged fluid, 
represented by the three-form field strength, and the black hole. If a neutral black 
hole is placed within a charged fluid, then some of the fluid flows into the black hole 
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Figure 1: Large r behavior of numerical black hole solutions. The figure plots the amplitude 
of the zero mode 70 at r = 10 6 as a function of h(ro), the electric field at the horizon, for 
a black hole with r$ = 10. The curve has a single zero at which the zero mode amplitude 
vanishes, and this uniquely determines a value of h(ro) for which the black hole is asymptotic 
to the Lifshitz fixed point geometry. 



and makes it charged. A static geometry describes an equilibrium configuration of 
the charged fluid outside a charged black hole, where the electric repulsion from the 
charge of the black hole precisely balances the gravitational pull. 

Let us assume that there is a non-degenerate horizon at r = r . The g tt component 
of the metric should then have a simple zero and the g„ component a simple pole at 
the horzion. If we further assume that the electric field h(r) has a finite value at the 
horizon, we find that the charged fluid density must go to zero at the horizon. This 
is in line with the equilibrium argument in the previous paragraph. 

With these assumptions, we can develop a near-horizon expansion of the various 
fields, 

f(r) = y/r - r (/ + fi(r - r ) + f 2 (r - r ) 2 + ...), 

g(r) = , l (g +gi(r - r ) + g 2 {r -r ) 2 + . . .) , 

j(r) = y/r - r (j + ji(r - r ) + j 2 (r - r ) 2 + ...), 

h(r) = h + h 1 (r-r ) + h 2 (r-r ) 2 + .... (29) 

Inserting this into the equations of motion and working order by order in r — r one 
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Figure 2: The charge of an asymptotically Lifshitz black hole is uniquely determined by the 
horizon area. The figure plots h(ro), the electric field strength at the horizon, as a function 
of ro, the area coordinate at the horizon. The dashed curve is the upper bound on ho in 
equation (34). The bound is saturated in the small black hole limit. 



obtains relations between the various constant coefficients, 

3/2 

9o = 7f=, , . f . (3°) 



v / FWi' 



Jo = , 2k0Vr ° , (31) 
2 h r ((h% - 3)r 2 - 1) 

~ — (5 - hiyi + 1 — ' { } 

/ ((6feg - 52feg + 100)r 4 + (45 - ll^)r 2 + 5) 
Jl 2r ((5 - hl)rl + l) 2 ' ^ 



We notice that for a given black hole size, 
field strength at the horizon, 



r , there is an upper bound on the electric 



(34) 




beyond which g and jo would be complex valued. 

We now use the expansion to generate initial values for the numerical integration 



of the equations of motion (17) - (20), starting close to the horizon and integrating 



outwards in r. We have a two-parameter family of initial data, using r and ho as 
the independent variables in the expansion. At first sight, fo also appears to be an 
independent free parameter but this is not really the case. Equation (17) is linear 
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Figure 3: A large black hole with ro = 20. The figure on the left shows the metric 
functions /(r) and g(r), while the figure on the right shows the electric field /(r) and the 
charge density j{r). 



in / so the overall normalization of f(r) is not determined. As discussed above, 
the solutions that are of most interest to us are those where / goes to a constant 
asymptotically, f(r) — > as r — > oo. Such a solution can always be normalized to 
/ — > 1 by dividing through by . We can therefore set fo = 1 in our numerical runs 
and take care of the normalization of / at the end of the day. 

The next step is to look systematically for solutions that describe black holes in an 
asymptotically Lifshitz background. A convenient way to conduct the search is to fix 
r and then prepare a sequence of initial data for different values of ho- For each set 
of initial values, the equations of motion are integrated from the near-horizon region 
out to sufficiently large r so that the non-zero eigenmodes have decayed away. One 
then looks for a zero mode of the form (24) and notes how the amplitude 70 varies 
as a function of ho, for a given value of ro- The result of this procedure is shown in 
Figure [TJ There is a unique value of h , for which the amplitude of the zero mode 
vanishes, and this corresponds to the charge of an asymptotically Lifshitz black hole 
with horizon at r = r . 

This can all be repeated for different sized black holes and Figure [2] shows the 
critical value of ho as a function of ro- This figure nicely summarizes the numerical 
results of this subsection. Starting from a two-parameter family of initial data, we 
have found a one-parameter family of black holes geometries that have the correct 
asymptotic behavior to be sitting in a Lifshitz background. The upper bound (34) on 
ho is given by the dashed curve in the figure. The bound is saturated, i.e. the black 
hole is becoming extremal, in the small black hole limit. 

The actual metric and gauge field configurations for two such black holes, a large 
one and a small one, are shown in Figures [3] and [4] respectively. 
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Figure 4: A small black hole with ro = 0.1. The figure on the left shows the metric 
functions f(r) and g(r), while the figure on the right shows the electric field /(r) and the 
charge density j(r). 

2.2 Lifshitz stars 

In this subsection we turn our attention to a different kind of localized object in 
an asymptotically Lifshitz background, one described by a smooth geometry with no 
horizon or curvature singularity anywhere in spacetime. It has a spherically symmetric 
concentration of the charged fluid with j(0) taking a finite value. Since the fluid 
density is finite at the origin the electric field strength must vanish there, which 
translates into h(0) = 0. If we also impose that the g u and g rr components of the 
metric be finite at r = we arrive at the following expansion at small r, 

j(r) = Jo (l-^(l + 2j 2 )r 2 + ^(103 + 152j 2 + 36j 4 )r 4 + •••), (35) 

h(r) = JorQ-^(6 + Jo 2 )r 2 + _±_( 5 28 + 96j 2 + 7j 4 )r 4 + . . . Y (36) 

g(r) = r(l + ^(-5 + j 2 )r 2 + ^(375 - 146j 2 -9j 4 )r 4 + ...^, (37) 

f(r) = ^(l + ^5 + 2j 2 )r 2 + ^(-125 + 16j 2 + 4j 4 )r 4 + ...^). (38) 

This can be used to generate a one-parameter family of initial value data, correspond- 
ing to different charge densities jo at the origin. We start the numerical integration at 
small r and integrate outwards to a large value of r, where the non-zero eigenmodes 
have decayed away. Figure [5] displays the amplitude of the remaining zero mode as a 
function of jo. It reveals a discrete set of 'magic' values of jo f° r which the zero-mode 
vanishes and the geometry is asymptotic to the Lifshitz fixed point geometry. 

We refer to the configurations with magic jo values as Lifshitz stars. They occur 
at j = 3.59, 21.8, 1.34 x 10 2 , 8.18 x 10 2 , 5.05 x 10 3 , 2.99 x 10 4 , . . .. These magic 
values can be seen in Figure [5] and the shape of the curve in the figure suggests that 
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Figure 5: Large r behavior of the non-singular solutions discussed in Section 2.2 The 
figure indicates the amplitude of the zero mode 70 at r = 10 4 as a function of jo, the charge 
density at the origin. The zeroes correspond to Lifshitz stars which are asymptotic to the 
Lifshitz fixed point geometry. We plot the product j'070 rather than 70 itself. This does not 
affect the location of the zeroes but makes the plot more readable. 



the sequence continues. We eventually run out of numerical precision when we go to 
higher j values. 

The metric and gauge fields of the Lifshitz star with the lowest magic value of jo 
are shown in Figure 16} 



3 Black hole thermodynamics 

In the original AdS/CFT correspondence, finite temperature is studied by replac- 
ing the AdSs part of the ten-dimensional spacetime by a five-dimensional AdS- 
Schwarzschild black hole [26] and we expect black holes to play a similar role here. 

The Hawking temperature can be obtained by going to a Euclidean metric and 
requiring regularity at the horizon, 

T H = ^rl (39) 
4?r g 



where we have used the near- horizon expansion (29 ) for the metric. The coefficients go 
and /o are easily determined from our numerical solutions for the metric and Figure [7] 
shows the Hawking temperature as a function of black hole size. 

We can immediately see important differences between the thermodynamic be- 
havior of asymptotically Lifshitz black holes and that of asymptotically AdS black 
holes. Let us first consider black holes that are small compared to the length scale 
L, set by the cosmological constant. Small AdS black holes behave much like ordi- 
nary Schwarzschild black holes in asymptotically flat spacetime. In particular, their 



10 



Figure 6: The Lifshitz star with the smallest allowed value of j(0), the charge density at 
the origin. The figure on the left shows the metric functions f(r) and g(r), while the figure 
on the right shows the electric field /(r) and the charge density j(r). 
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Figure 7: Numerical results for Hawking temperature as a function of black hole size. 
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Figure 8: Black hole entropy as a function of Hawking temperature. 



Hawking temperature increases as they get smaller leading to a thermodynamic in- 
stability. In the asymptotically Lifshitz case, on the other hand, the temperature is a 
monotonic function of black hole size and is rapidly falling at the smallest black hole 
sizes that our numerical calculations can handle. This supports our earlier claim that 
these black holes become extremal in the limit of vanishing black hole size. It also 
means that there is no analog of the Hawking-Page transition in the Lifshitz system. 
Black holes that are large compared to L satisfy simple scaling relations. The 



form of the near-horizon expansion (29) suggests that go 



,1/2 



and f 



-1/2 



for 



black holes with tq 3> 1, which in turn gives 

S oc T H , 



(40) 



where S = irrl is the Bekenstein-Hawking entropy (the corresponding behavior for 
large 3+1 dimensional AdS black holes is S oc T%). This is confirmed by our numerical 



results. We find that go ~ 0.57 r. 



1/2 



and fo ~ 1.96 r 
S ss 11.4 T H 



-1/2 



giving 



(41) 



for large asymptotically Lifshitz black holes. Figure [8] plots the black hole entropy 
as a function of temperature over the range of black hole sizes for which we have 
obtained numerical solutions. 



4 Wilson loops 

The bulk theory we are studying is conjectured to be dual to a boundary theory at 
the Lifshitz point. Following j20] we observe that the term K (V 2 0) in the Lifshitz 
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action (2J) can be written K | V x E\ , where Ei = eijdj(p automatically solves Gauss' 
law. Note also that the field is dimensionless under Lifshitz scaling (fll), K is a 
dimensionless coupling constant, and an additional term in the action of the form \E\ 
would be accompanied by a coupling with dimensions of mass squared. We conclude 
that the boundary theory can be viewed as a gauge theory in 2 + 1 dimensions with a 
dimensionless coupling constant and an unusual action. The dimensionless coupling 
suggests that the theory perhaps has some nice features in common with conventional 
gauge theory in 3 + 1 dimensions. 

With the dual gauge theory in mind, we introduce Wilson loops on the gravity side. 
The Wilson loops contain information about the force acting between 'quarks', i.e. 
particles that are charged under the gauge fields. The recipe given in [271 12H] involves 
hanging a string from the boundary, with the end points of the string representing 
the quarks. 

The action of the string for a rectangular Wilson loop, with initial and final 
Euclidean time separated by A, is given by 
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dx\\ f 2 g 2 r 2z ~ 2 ( ^) +pr 2z + 2 , 



(42) 



where we have put a = x and r = t in static gauge. Extremizing the action leads to 

j-2 r 2z+2 



z+l 
mm' min > 



(43) 



P2„2z+2 



where r min is the r coordinate of the midpoint of the hanging string and / m ; n = /(r min ). 
The boundary distance between the end points of the string is then given by 

dr 
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(44) 
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The energy of the string configuration is 
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(45) 



where we have regularized the expression by subtracting the contribution of two 
straight strings going from the boundary down to the horizon of the black hole at 
r = r . 
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Figure 9: Screening behavior at finite temperature. The figure on the left shows the 
boundary distance between the endpoints of a hanging string as a function of r m ; n , while 
the figure on the right plots the potential energy as a function of r m i n . 

At vanishing temperature there is no black hole and the above expressions simplify 
considerably. In this case /, g — > 1 and r — > 0, and the boundary distance between 
the endpoints of the strings is given by 



dr 1 2 f°°dy 1 _ 2v^r(|^) 



while the potential energy reduces to 

v - 1 1 A 1 

z 



(46) 




1 

a>L \ i J ^z Ir(^) 



(47) 



Here is a dimensionless coupling constant that we assume to be large. The energy 
V has dimensions of inverse length, hence the factor of 1/L, and the dependence 
on I can be traced to the unconventional scaling properties of the Lifshitz system. 
At z = 1 we recover the results of [271 EH]- At z = 2 it is tempting to relate the 
dimensionless coupling to the Lifshitz coupling K. One has to keep in mind, though, 
that our calculation can at best be expected to make sense at strong coupling whereas 
K is denned in a free theory. 

At non-zero temperature we have to resort to a numerical evaluation of the inte- 
grals, using our numerically evalutated metric as input. Figure [9] shows the boundary 
distance between the string endpoints as a function of r m i n and the potential energy 
of the hanging string, also as function of r min , for z — 2. The black hole is taken 
to be large so that I <C r . The figure shows qualitatively the same behavior as 
was found for 3 + 1 dimensional gauge theory in [291 E0]- At small separation the 
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Figure 10: The figure on the left shows the potential energy as a function of the boundary 
distance between quarks represented by the endpoints of the hanging string. The figure 
on the right plots the critical distance, £ c , where screening sets in, vs. temperature. The 
log- log graph shows a clear T~ 1 / 2 dependence. 

inter-quark potential is similar to the zero temperature potential but when the sepa- 
ration between the quarks becomes sufficiently large the gauge interaction is screened 
and the potential energy vanishes. The crossover occurs where the potential energy 
V becomes positive in Figure [9| signaling that at large separation the configuration 
minimizing the energy is simply two straight strings stretching from the boundary 
down to the horizon. Since there is no interaction between such strings the potential 
simply vanishes in this case. Figure [TUIplots the resulting potential as a function of the 
boundary distance between the string endpoints, and also shows how the crossover 
distance, where screening sets in, depends on temperature. The Lifshitz scaling is 
apparent in the £ c oc T" 1 / 2 falioff as opposed to the T~ l behavior seen in 3+1 dimen- 
sional gauge theory. 

We do not know whether the finite temperature behavior exhibited by the Wilson 
loops is in some way related to the ultra- locality discussed in [TBI EI] ■ We hope to 
return to this issue in future work. 

5 Conclusions 

In this paper we have explored further the recently proposed gravity dual description 
of Lifshitz type fixed points. We have mainly focused on the gravity side of the 
duality, finding non-trivial spacetime geometries that are asymptotic to the Lifshitz 
fixed point geometry, including black holes that provide a dual description of a Lifshitz 
system at finite temperature. 

It is by no means obvious how to incorporate the bulk metrics that we have 
found into a solution of ten dimensional string theory. It is nevertheless tempting 
to proceed under the assumption that such a construction can be found, and that a 
duality analogous to AdS/CFT actually exists. Alternatively, we can view the gravity 
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dual as a purely phenomenological description of the 2 + 1 dimensional physics. Either 
way, one is motivated to study the gravitational theory in more detail. 

Our numerical black hole solutions in principle contain all the information that is 
needed to calculate finite temperature correlation functions in the dual system but 
this requires rather delicate numerical analysis which we leave for future work. It 
also requires a better understanding of holographic renormalization for non-Lorentz 
invariant field theories [32J. 

It would also be interesting to make contact with recent work [33] on non-relativistic, 
non-abelian gauge theories which exhibit z = 2 quantum critical behavior. A gravity 
dual of a large-N limit of such a theory would be on firmer theoretical ground than 
the models studied here. 
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